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ABSTRACT 



X 

5_j ' Aims. To study the possibility of appearance of accelerated universe in scalar tensor cosmological models 



cd 



Methods. We consider scalar tensor theories of gravity assuming that the scalar field is non minimally 
coupled with gravity. We use this theory to study evolution of a flat homogeneous and isotropic universe. 
In this case the dynamical equations can be derived form a point like Lagrangian. We study the general 
properties of dynamics of this system and show that for a wide range of initial conditions such models 
lead in a natural way to an accelerated phase of expansion of the universe. Assuming that the point like 
Lagrangian admits a Noether symmetry we are able to explicitly solve the dynamical equations. We study 
one particular model and show that its predictions are compatible with observational data, namely the 
publicly available data on type la supernovae, the parameters of large scale structure determined by the 
2-degree Field Galaxy Redshift Survey (2dFGRS), the measurements of cosmological distances with the 
Sunyaev-Zel'dovich effect and the rate of growth of density perturbations 

Results. It turns out that this model have a very interesting feature of producing in a natural way an epoch 
of accelerated expansion. With an appropriate choice of parameters our model is fully compatible with 
several observed characteristics of the universe 



Key words, cosmology: theory - cosmology: dark energy - quintessence - large-scale structure of Universe- 
scalar 
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1. Introduction 

Recent observations of the type la supemovae, Gamma Ray Bursts, and CMB anisotropy indicate 
that the total matter-energy density of the universe is now dominated by some kind of dark energy 
JRiess & al. 1998llRiess 2000IIRiess & al. 2004^ . The origin and nature of this dark energy is not 
yet known d^eldovich 1^67l|Wemberg m')\ . 

In the last several years a new class of cosmological models has been proposed. In these 
models the standard cosmological constant A-term is replaced by a dynamical, time-dependent 
component - quintessence or dark energy - that is added to baryons, cold dark matter (CDM), 
photons and neutrinos. The equation of state of the dark energy is given by w^ = p^/p,/,, where 
p^ and p0 are, respectively, the pressure and energy density, and - 1 < w^ < 0, what implies 
a negative contribution to the total pressure of the cosmic fluid. When w^ - -1, we recover 
a constant A-term. One of the possible physical realizations of quintessence is a cosmic scalar 
field ( ICaldwell, Dave & Steinhardt 1998t , which induces dynamically a repulsive gravitational 
force, causing an accelerated expansion of the Universe, as recently discovered by observations 
of distant type la supernovae (SNIa) JRiess & al. 19981 IRiess & al. 2004t and confirmed by the 
WMAP observations ( |Spergel & al. 2003> . 

The existence of a considerable amount of dark energy leads to at least two theoretical 
problems: 1) why only recently dark energy started to dominate over matter, and 2) why dur- 
ing the radiation epoch the density of dark energy is vanishingly small in comparison with 
the energy density of radiation and matter (fine tuning problem). The fine tuning problem can 
be alleviated by considering models of dark energy that admit so called tracking behavior 
l |Steinhardt, Wang & Zlatev) . In such models, for a wide class of initial conditions, equation 
of state of dark energy tracks the equation of state of the background matter and radiation 
jSteinhardt, Wang & Zlatev] [Zlatev, Wang & Steinhardt 1999| . All these circumstances stimu- 
lated a renewed interest in the generalized gravity theories, and prompted consideration of a 
variable A term in more general classes of theories, such as the scalar tensor theories of grav- 
ity ( | Perrotta, Baccigalupi & Matarrese 2000) . One of the additional advantages of appealing to 
these theories is that they open new prospective in the scenario of a decaying dark energy, since 
the same field that is causing the time (and space) variation of the dark energy is also causing the 
Newton's constant to vary. 

In this paper we consider cosmological models in a non minimally coupled scalar tensor 
theory of gravitation. We assume that the universe is homogeneous and isotropic and its geometry 
is described by the Friedman-Robertson- Walker line element. For the reason of simplicity we 
consider only the case of flat universe. In this case the scalar field depends only on time and the 
dynamical equations that describe evolution of the geometry and the scalar field can be derived 
form a point like Lagrangian. We derive the general set of dynamical equations and discuss their 
basic properties. When we require that the point like Lagrangian admits an additional Noether 
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symmetry the dynamical equations can be explicitly integrated. In particular we consider a model 
with the scalar field potential of the form V((f>) = Vo^'* and we analyze its dynamics. Finally to 
compare predictions of our model with observations we concentrate on the following data: the 
publicly available data on type la supernovae, the parameters of large scale structure determined 
by the 2-degree Field Galaxy Redshift Survey (2dFGRS), and measurements of cosmological 
distances with the Sunyaev-Zel'dovich effect. We show that our model is compatible with these 
observational data. 

Our paper is organized as follows: in Sec. 2 we present our model and discuss its basic 
properties. In Sec. 3 we confront predictions of our model with observational data. Sec. 4 is 
devoted to the discussion of evolution of density perturbations in our model and finally in Sec. 5 
we present our conclusions. 

2. Model description 

Let us consider the general action of a scalar field non minimally coupled with gravity when 
there is no coupling between matter and cp 

Ji^ ^ ^Fg (f{^)R + ^/>,;,^,v - V{(P) + £j d'x, (1) 

where F(ip), V((p) are two generic functions representing the coupling of the scalar field with 
geometry and its potential energy density respectively, R is the curvature scalar, —g^^(p^^(py is 
the kinetic energy of the scalar field (p and £.,„ describes the standard matter content. In units 

such that ^jiGn - fi - c - 1, where Gn is the Newtonian constant, we recover the standard 

1 , 

gravity when F - — , while in general the effective gravitational coupling Geff — —jp- Here we 

would like to study the simple case of a homogeneous and isotropic universe, what implies that 

the scalar field depends only on time. It turns out that in the flat Friedman-Robertson- Walker 

cosmologies, the action in Eq. Q reduces to the pointlike Lagrangian^ 

£ = 6Faa^ + 6F'4>a^d + a^ I -(^^ - ViijM - Da"^<>'"" , (2) 

where a is the scale factor and prime denotes derivative with respect to 0, while dot denotes 
derivative with respect to time. Moreover, the constant D is defined in such a way that the matter 
density p,„ is expressed as p,„ - D^aojay'^, where 1 < y < 2. The effective pressure and energy 
density of the 0-field are given by 

P^ = ^0' - V{4>) - 2(F + 2HF) , (3) 

P^^^^^ + V(<f>) + 6HF, (4) 



' We use the expression pointlike to stress that the field Lagrangian obtained from Eq. Q can be con- 
sidered as defined in the minisuperspace where the remaining two variables (a,(^) depend only on the 
cosmological time t and so they can be considered as describing a mechanical system with two-degrees of 
freedom. 
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a 
where // = - is the Hubble constant. These two expressions, even if not pertaining to a conserved 

a 

energy-momentum tensor, define an effective equation of state Wa, - — , which drives the late 
time behavior of the model. The field equations derived from Eq. (|5} are the same ones which 
would come from the field equations derived from Eq. Q when homogeneity and isotropy is 
imposed, that is 



2F I 3 3 / 



(5) 



2H + 3//2 = -L(p^ + p^) , (6) 

and 

if) + 3H^ + 6{H + 2H^)F' + V = . (7) 

The scalar tensor theory with the scalar field non minimally coupled to gravity provides a wide 
framework to study cosmological models. The function Fi^i) that describes the coupling between 
the scalar filed and gravity is influencing not only the evolution of the cosmological scale factor 
and the scalar field itself but also determines the strength of gravitational interactions. In this 
paper we consider only the simple case of a homogeneous and isotropic flat universe filled in with 
a scalar field (quintessence) and pressureless matter, i.e., p„ - (dust). That is, strictly speaking, 
our model describes the evolution of the universe only after the matter-radiation decoupling. 
The scalar field which appears in our model is treated as quintessence. To discuss the changing 

influence of the scalar field on the evolution of the universe and the effective equation of state of 

(jp- F 

quintessence, we set x = — , and xi = — , we find that 

p , x-\+ Ixx (^ -I- 2//) - 2ii 

W0 = — = ■ — 777 • ^^) 

P0 x+\ + btix\ 

The parameter x measures the ratio of the kinetic energy relative to the potential energy of the 
scalar field. In the non minimally coupled case it is possible to invert the Eq.ljSj and we get that 

P^(1+V + 2(f-//f) 

X - — ^^ (9) 

P0(1-W0)-2(f + 5//f) 

1 + W0 

When the coupling function F is a constant we recover the relation x - that holds in 

1 - W0 

the minimally coupled case. When w^ < the quintessence contributes negative pressure, this 
occurs when 



x<\+2 



xi-xi\—+lH 



(10) 



while the inequalities 

x>xx-x,\^-+5H\, (11) 

X < xi-xA— +5H\ , (12) 
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correspond to the standard quintessence (w^ > -1) and superquintessence respectively (w^ < 
-1). If both F((f)) and V{(p) are known, it is possible to describe, in the parameter space, the 
transition between standard quintessence and superquintessence. Let us now introduce the con- 
cept of an effective cosmological constant A^,//. Using Eq.(|5|i it is natural to define the effective 

cosmological constant as Ae/f - -7— . With this definition we can rewrite Eq.(|5ji as 

zF 

3H^ = G,ffp„,+Kff. (13) 

Introducing the standard Omega parameters by 

o P"- o _ ^"ff _ P-A 

'"" 6F//2' *" 3//2 " 6F//2' 

we get that as usual 

Q„, + 0^ = 1. (14) 

^From the definition of p^ and p^ and the generalized Klein-Gordon equation it follows that 

P^ + 3H[p^+P^)^-6H^F, (15) 

and 

Aeff + GeffPm = -3HGeff [p^ + p^) ■ (16) 

These equations play a role of the continuity equations forp^ and Af//. 

We will show later on that asymptotically for large time Geff tends to zero, then w^ determines 

the late time scaling of A^//, and actually we have that 

^^-37?(l+wJ, (17) 

where H is the asymptotic value of H. It turns out that if at large f, H is constant and + 0, 
than asymptotically w^ — > -1. Let us note that such a transition to the asymptotical value is 
responsible for the accelerated expansion. Even if we start from a dust or other stiff equation of 
state W0 converges towards -1, as it is shown in Fig^ Actually it turns out that together with 
the superquintessence region w^ < -1, (SQ), there exists also a superquintessence connected 
region (SQC) with an equation of state w^ > 0. The blue straight line corresponds, indeed, to 
the values of parameters such that we observe neither a superquintessence expansion, nor a stiff 
matter behaviour. 

2.1. Exact solutions through Noether theorem 

To find exact solutions in the framework of the non minimally coupled models we assume that the 
pointlike Lagrangian possesses a Noether symmetry (for a detailed exposition of this technique 
which we closely follow here and interesting examples see ( |Capozziello & al. 1996| ). This means 
that we require the existence of a Noether vector field along which the Lie derivative of the 
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Fig. 1. Rate of change of the equation of state as measured by w^ versus the w^ parameter. We 
see that together with the superquintessence region w^ < -I, (SQ), there appears also a su- 
perquintessence connected region (SQC) with an equation of state w^ > 0. The blue straight line 
corresponds, indeed, to values of parameters such that we observe neither a superquintessence 
expansion, nor a stiff matter behaviour 

Lagrangian is zero. This requirement restricts the possible coupling and the form of the potential. 
In fact the additional Noether symmetry exists when 



V^Vo(F(<p))'"-'\ 

where Vo is a constant and 

,, 3(5+1) 
p(s) 



(18) 



(19) 



2s + 3 
where i is a real number, and when the coupling F((p) satisfies the following differential equation 



diF"F^ + diF"^ + d-iF'^F + d^F^ = 0, 



and the coefficients are functions of the parameter s, 

2s + 3 



di ^3s{s+ l){s + 2), 

d^ ^--(s+l)(Ss^ + l6s + 3), 



d^ 



s(2s + 3)^ 
12 



A particular solution which indeed turns out to be quite interesting is of the form 

F = ^(i)(0 + <p^f , 
where 



(2s + 3f 



48(i + l)(s + 2) 



(20) 

(21) 

(22) 
(23) 

(24) 
(25) 
(26) 
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and 00 is a constant. The parameter s labels then the class of Lagrangians which admit a Noether 
symmetry. Let us note that the form of the coupling given by ( I25> is quite relevant from the 
point of view of fundamental physics. Not all the values of s are allowed however From the 
expressions (I18> and ( I26> it follows that the cases when s - -I, s = -3/2, s - -2 are 
special and they should be treated independently. Through Eq. ( I18> it is possible to recover 
the inverse power-law potentials, while the case .s = as we will see later is special and it 
corresponds to the square hyperbolic sine potential. Both these potentials are usually assumed 
ad hoc to obtain certain asymptotic behavior of the energy density of the quintessence field 
dPeebles & Ratra 1988l|Urena-Lopez & Matos 2000| l, while here they emerge naturally from the 
imposed Noether symmetry. Once the general solution of Ea. J20> is found, what automatically 
specifies the form of the potential, and a generic value of s is considered, it is possible to ex- 
plicitly find the Noether symmetry and to introduce new dynamical variables associated with 
this symmetry (for details see ( |Capozziello & al. 1996| l). Using the new dynamical variables it 
is then possible to solve the corresponding Lagrange equations and finally by inverting them we 
obtain the sought after a(t) and (f>(t). The final result can be written in the form 

a(t) = A(s){B(sy- ^f- + — f T^ , (27) 

\ 6(^-1-6) Eo/ 

Vq 3 D\ ^' (2«3)2 

cpit)^Cis)\--^B(s)t- + —\ t-^^-4,0, (28) 

\ 7(s) So/ 

where D is the matter density constant, Zq is a constant of the motion resulting from the Noether 
symmetry, Vq is the constant that determines the scale of the potential, (po is a constant that 
determines the initial value of the scalar field and the other constants A{s), B(s), C(s), and y( s) 
are given in the appendix. As it is apparent from Eq.j27> and Eq.j28> for a generic value of s 
both the scale factor a(t) and the scalar field (p(t) have a power law dependence on time. It is also 
clear that there are two additional particular values of s, namely s - and s = -3 which should 
be treated independently. In this paper we concentrate on the case s - -3 which provides, as will 
be shown shortly, an interesting class of models. 

2.1 .1 . The case of quartic potentials: analysis of the solution 

When s = -3 the general solutions given by J27> and J28> lose their meaning. To find a(t) and (p{t) 
in this case it is necessary to use the general procedure as described in ( |Capozziello & al. 1996^ . 
From Eg.ltrSt and Eq.(|53) it follows that in this case F - ^<^^, and V(<f>) = Vo<f>'^ where we have 
set (/>o = and Vo denotes a constant. This case is particularly interesting since the resulting self- 
interaction potential is used in finite temperature field theory. In fact, it seems that the quartic 
form of the potential is required in order to implement the symmetry restoration in several Grand 
Unified Theories. Just for this reason we limit our analysis to this special case, and we will show 
that it provides an accelerated expansion of the universe. We reserve to a forthcoming paper the 
study of the other cases. As in the general case once we have the functions F{<p) and y(0) that 
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allow a Noether symmetry it is possible to explicitly find it and effectively solve the Lagrange 
equations. Performing this procedure we finally get 

ait) = Q-oe^ [(e"' ' - l) + ai f + as] ' , (29) 

m = 00 J , „, '"'' ^^ , (30) 

y (e"^ - 1) + 0-2 f + as 

where ao, 01,020:3 and ^o are integration constants. They are related to the initial matter density, 
D, and the scale of the potential Vq hy D - j^a^cfy^^a 102, which implies that they cannot be zero. 
The case Vq - Q has to be treated separately. It turns out that the constants a^, 0o and oq have 
an immediate physical interpretation: oq, and 03 are connected to the value of the scale factor at 
f = 0, actually fl(0) - (03)^01). Moreover 03 can be selected in such a way that at a sufficiently 
early epoch the universe is matter dominated. This requires that 03 be sufficiently small, for 
example, that 03 e [0.001,0.01]. The constants 0-3 and 0o are connected to the initial value 

00 

of the scalar field cf>(0) - —= and therefore they determine the initial value of the effective 

\a3 



16a3 
gravitational constant Gf//(0) = r-. Let us note that an attractive gravity is recovered when 

00 is a pure imaginary number Without compromising the general nature of the problem we can 

set, for example (po - i. This choice does not violate the positivity of energy density of the scalar 

field or the weak energy condition. To determine the integration constants o-i, 02, 0:3 and 0o we 

follow the procedure used in JDemianski & al. 2005> . and we set the present time to - I. That is 

to say that we are using the age of the universe, fo, as a unit of time. Because of our choice of time 

unit the expansion rate H(t) is dimensionless, so that our Hubble constant is not the same as the 

Hq that appears in the standard FRW model, measured in kms^^Mpc^^: we then set Hq - //(I). 

Using (29) we get 

H(l) = //o = --^ + 



3 3 e"' + a2 + a3 - 1 ' 

which we use to find q'2 in the form 

e"'(ai - 3//0) + (ai - l)(3//o + ai) 

Q'2 = — . 

3//0 + ai - 2 

With this choice of time the scale factor, the scalar field and the expansion rate assume the final 

form 

-T I .«.' ^ g"''(«i-3J/o) + (a3-l)(3J?o + gi) ^ ^ ^^ _ ^^'^' 



(31) 



3//0 + Cl 
(pit) = 00 _ 



\ gff,/ _ e"! (a-i -3Ho)+(g3-l)(3Ho+ffi ) ^ ^ _ -^ ' 

\ 3Ho+ai-2 

Hit) = [ay (ai + 2e"') + e"''ai [3% + ai - 2) + 3Hq (a, - 2e^ + 2) (32) 

- {a\ + 3Ho (o-i + 2)) 03 + [(-1 - 4) a\ + 030] + 3 (-1 + e" + as) //oo-i] t] X 
{{3e"" (3//0 + ai - 2) - {03 - 1) (3//o(f - 1) + o-^f - 1) + 2) + 6^(0-1 - 3//o)f}"' . 
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Fig. 2. Plot of log,o<5G,„ = logi, 
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Fig. 3. Time dependence of w^, for two different values of the parameters. Actually the solid curve corre- 
sponds to cfi = 2.5, and Hq = 0.95, while the dashed one corresponds to a^ = 3, and Hq = I. We see that 
even though they both produce accelerated expansion, only the second one gives rise to super acceleration. 



Let us remind that t is now varying from to 1 and t - 1 corresponds to the present moment. 
The parameters Hq and ai admit a simple physical interpretation. Actually Hq is the present 
value of the Hubble constant measured in our unit of time, while o-i drives the early time and the 
asymptotical behavior of a{t) and (^(f). For f <k — we have 



a(t) 



(a\ + a2 — )t + as 



(p(t) ~ [{ai + a2 - aia3)t + a^Y 



(33) 
(34) 



Later at larger t,a{t) reaches an intermediate stage, when it evolves as ait) ~ Pe^~ (dumped 
dust) , and has a de Sitter behavior a(t) ~ e~ for f — > co. 

It is interesting to note that^ w^ is representing an equation of state, in the usual sense, of 
the effective cosmological constant Af// . In Figs. (l3l4l5> we show the time dependence of w^: 
we see that this equation of state can admit a superquintessence behavior (w < -1), just as an 
effect of the transition toward w — > - 1 . We note that in the remote past as in the far future w^ is 



^ Moreover, we will show in the following that in order to fit the observational data H = ^ has to be of 
the same order of magnitude as Hq, i.e. roughly ai « 3Hq. This implies that this model even if formally 
depending on two parameters is, as a matter of fact, very sturdy and depends mainly on the Hubble constant. 
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Fig. 4. Parametric plot of w^ as a function of the acceleration; we see the transition from the accelerated to 
the decelerated expansion. 
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Fig. 5. Parametric plot of w^ with respect to the deceleration parameter q = —^\ again we see the transition 
from the accelerated to the decelerated expansion. Moreover, we see that if we consider the value q ^ -0.68 
obtained from the SNIa Gold Sample IRiess & al. 20041 we obtain w^ < -0.65 



constant: it mimics an almost dust equation of state (w^ ~ 0) in the far past and asymptotically 
behaves as a bare cosmological constant (w^ — > -1) as f — » oo. Let us also mention that since 
both p^ and p^ depend on F{^) through its time derivative and asymptotically (pit) ~ constant we 
asymptotically recover the minimally coupled case, with 






(35) 
(36) 



However, before reaching this asymptotic regime the total energy density p^ is dominated by the 
coupling term 6HF. In conclusion of this section we present the traditional plot logp^ - log a 
compared with the matter density (see Fig. (|6|l). Interestingly we see that p^ tracks the matter 
during the matter dominated era ( actually a(t) oc ff , and w^ ~ 0), and becomes dominant at 
late time. In Fig.Q we plot the redshift behavior of the quintessence sound velocity; we see that 
Cj < 1. 
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Fig. 6. Plot of log,QP(j versus logjQfl. The vertical bar marks logigao. The solid blue straight line indicates 
the log-log plot of p„ versus a. 
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Fig. 7. The redshift behavior of the quintessence sound velocity Cj 



15 


\ 




glO 


\x 




^ 5 


^N. 







\^\^ 




^~^ 






^\^ 



-4 -3 -2 -1 
Log (a) 

Fig. 8. Plot of logiQp^ versus logjoO when there is asymptotic freedom at f ^ . The vertical line marks 
logjoflo- The solid blue straight line indicates the log-log plot of p^ versus a. 



2.1 .2. A special case: asymptotic freedom at ? ^ 

In this section we consider a special case when lim,^o Ge/f - 0, that is when lim,^o = oo we 
have a sort of asymptotic freedom at f = 0. First of all we note that such a case can be reached 
by setting 0-3 = 0, which also implies that a(t = 0) = 0, and the expressions for a{t) and <p(t) 
become simpler. In comparison with the case as ^ the main difference concerns the behavior 
of the density p^ with respect to the matter density, as shown in Fig. (|8|i. Actually we note that 
the coupling F((p) diverges as f — > 0, with its derivatives, more rapidly than p,„. However p,„ 
always dominates over the scalar field contribution to the density, that is over p^. It should be 
remembered that the early epoch f « does not belong to the physical time domain of our 
model, since we are neglecting the contribution of radiation, and therefore the time behavior of 
Geff during the intermediate period is unknown. This is also the reason why we are not using 
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Fig. 9. The same as in Figs <4l5t when there is asymptotic freedom at / ^ . 



any constraint on Geff from the nucleosynthesis. With respect to the other characteristic features 
this case is similar to that discussed above, as shown in the Fig. (|5J 

3. Observational data and predictions of our models 

Above we discussed some general properties of our scalar field model of quintessence, stressing 
how it provides a natural mechanism for the observed accelerated expansion of the universe. To 
test viability of our model we compare its predictions with the available observational data. We 
concentrate mainly on two different kinds of observational data: some of them are all based on 
distance measurements, as the publicly available data on type la supernovae, the measurements 
of cosmological distances with the Sunyaev-Zel'dovich effect and radio-galaxies data, other are 
instead connected with the large scale structure, such as the parameters of large scale structure 
determined by the 2-degree Field Galaxy Redshift Survey (2dFGRS), and the gas fraction in 
clusters. Before we begin our analysis we note in Fig.[21that for our model the transition redshift 
from a decelerating to an accelerating phase in the evolution of the universe falls very close to 
z - 0.5, in agreement with recent results coming from the SNla observations (IRiess & al. 2004> . 



3.1. Constraints from recent SNla observations 

In recent years the confidence in type la supernovae as standard candles has been steadily grow- 
ing. Actually it was just the SNla observations that gave the first strong indication of an acceler- 
ating expansion of the universe, which can be explained by assuming the existence of some kind 
of dark energy or nonzero cosmological constant JSchmidt & al. 1998t . Since 1995 two teams 
of astronomers have been discovering type la supernovae at high redshifts. First results of both 
teams were pubhshed by Schmidt & al. (1998) and Perlmutter & al. (1999). Recently the High-Z 
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Fig. 10. Observational data of the Gold Sample of 
SNIa (Riess et al. 2004) and the SNLS dataset 
<Astier & al. 20051 fitted to our model. The solid curve 
is the best fit curve with Hq 
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which corresponds to Hk^^ = 0.^4^^'^^. We also get 



/z = 0.68+;j^. 




Fig. 11. The same as in the Eq. <10> , but zooming on the high 
redshift SNIa. 




Fig. 12. Behavior of the second derivative of the scale factor. Please note that the transition from a decel- 
erating to an accelerating expansion occurs close to z = 0.5, as predicted by recent observations of SNIa 
z, = 0.46 ± 0.13 JRiess & al. 20041 . 



SN Search Team reported discovery of 8 new supernovae in the redshift interval 0.3 < z < 1.2 and 
they compiled data on 230 previously discovered type la supernovae l |Tom"y et al. 2001^ . Later 
Barris & al. (2004) announced the discovery of twenty-three high-redshift supernovae spanning 
the range of z = 0.34 - 1.03, including 15 SNIa at z > 0.7 . 

More recently Riess & al. (2004) announced the discovery of 16 type la supernovae with the 
Hubble Space Telescope. This new sample includes 6 of the 7 most distant (z > 1.25) type la 
supernovae. They determined the luminosity distance to these supernovae and to 170 previously 
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reported ones using the same set of algorithms, obtaining in this way a uniform "Gold Sample" 
of type la supernovae containing 157 objects. Finally just recently the Supernova Legacy Survey 
(SNLS) team presented the data collected during the first year of the SNLS program. It consists 
of 71 high redshift supernovae in the redshift range z e [0.2, 1]. This new SNLS sample is char- 
acterized by precise distance measurements of aU the 7 1 supernovae, so it can be used to build 
the Hubble diagram extending to z = 1 . The purpose of this section is to test our scalar field 
quintessence model by using the best SNla data sets presently available. As a starting point we 
consider the gold sample compiled by (IRiess & al. 2004> to which we add the SNLS dataset. To 
constrain our model we compare through a x^ analysis the redshift dependence of the observa- 
tional estimates of the distance modulus, ji - m - M, to the corresponding theoretical values. 
The distance modulus is generally defined by 

m-M = 51ogDz.(z) + 51og(-^) + 25, (37) 

where Ho is the standard Hubble constant, measured in kms^^Mpc^^, m is the appropriately 
corrected apparent magnitude including reddening, K correction etc., M is the corresponding ab- 
solute magnitude, and Di is the luminosity distance in Mpc. However, in scalar tensor theories 
of gravity it is important also to include in the Eq. (j^) corrections, which describe the efifect of 
the time variation of the effective gravitational constant Gefj on the luminosity of high redshift 
supernovae. Actually, if the local value of G^ff at the space time position of the most distant 
supernovae differs from On, this could in principle induce a change in the Chandrasekhar mass 
Mch cc G^i. Some analytical models of the supernovae light curves predict that the peak lumi- 
nosity is proportional to the mass of nickel produced during the explosion, which is a fraction of 
the Chandrasekhar mass. The actual fraction varies in diff'erent scenarios, but always the physical 
mechanism of type la supernovae explosion relates the energy yield to the Chadrasehkar mass. 
Assuming that the same mechanism for the ignition and the propagation of the burning front is 
valid for SNIa at high and low redshifts, it turns out that the predicted apparent magnitude will 
be fainter by a quantity jGaztaiiaga & al. 2002 ) 

AMG = ^log(^). (38) 

4 \GeffJ 

Taking this into account the distance modulus becomes 

m - M = 5 logDiiz) + 5 log(— ) -H 25 -H AMc- (39) 

The presence of this correction actually allows one to test the scalar tensor theories of gravity 
l |Gaztanaga & al. 2002 | lUzan 20(Bt using the SNIa data. For a general flat and homogeneous 
cosmological model the luminosity distance can be expressed as an integral of the Hubble func- 
tion as follows: 

Dl(z) = -^(1+z) [ T^,d^' (40) 

no Jo /^(O 

where H(z) is the Hubble function expressed in terms of z = ao/a(t) - 1. Using Eqs. J38l l and 
( I40t . which in our case can be integrated only numerically, we construct the distance modulus 
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and perform the x^ analysis on the complete data set. We obtain x%j =1.7 for 230 data points, 
and the best fit value is ^o = l-0!oo4, ^^i = 2.9+°-^, which coiTesponds to Qa^^ = ^■'^^-om- We 
also get h = 0.68^ qj. Moreover we checked that any value of a^ e [0.001, 0.01] does not aff'ect 
the determination of distances, so that in the following we set 03 = 0.001, being confident that 
such a choice does not alter the main results. In Fig JlOt we compare the best fit curve with the 
observational data sets. 

3.1 .1 . Dimensionless coordinate distance test 

After having explored the Hubble diagram of SNIa, that is the plot of the distance modulus as 
a function of the redshift z, we want here to follow a very similar, but more general approach, 
considering as cosmological observable the dimensionless coordinate distance defined as : 



f- 



y(^)- I j^/(- (41) 

It is worth noting that y(z) does not depend explicitly on h so that any choice for h does not alter 
the main result. Daly & Djorgovski ( |Daly & Djorgovski 2004| l have determined y(z) for the SNIa 
in the Gold Sample of Riess et al. tRiess & al. 20041 1 which represents the most homogeneous 
SNIa sample available today. Since SNIa allows to estimate D^ rather than y, a value of h has to 
be set. Fitting the Hubble law to a large set of low redshift (z < 0.1) SNIa, Daly & Djorgovski 
l |Daly & Djorgovski 2004| i have found that : 

h = 0.66 + 0.08 km s"' Mpc"' , 

which is consistent with our fitted value h = 0.66^ q^. To enlarge the sample, Daly & Djorgovski 
added 20 further points on the y(z) diagram using a technique of distance determination based 
on the angular dimension of radiogalaxies ( |Daly & Djorgovski 2004^ . This extended sample that 
spans the redshift range (0.1,1.8) has been obtained by homogenizing dififerent kinds of measure- 
ments, affected by different systematics, so that the full sample may be used without introducing 
spurious features in the y{z) diagram. However before using the dimensionless coordinate dis- 
tance test we do not use the Daly & Djorgovski database directly, but, for the supernovae Gold 
Sample we first converted the distance modulus into yizd using the Eq. d^ and the relation 

DL = -^(l+z)y(z). 
no 

Here again Hq is the standard FRW Hubble constant. To determine the best fit parameters, we 
define the following merit function : 

vC?:: ITi _ Hn^ — v.- 

(42) 



x\auHo)^^Yj 



N-3^__^ 



y(zr,ai,Ho)-yi 



We obtain ;k'^^^ = 1.19 for 186 data points, and the best fit value is Hq - 0.98^qq-3, qti = 2.5^q2. 
In Fig ( I13> we compare the best fit curve with the observational data set. Daly & Djorgovski 
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Fig. 13. Observational Daly & Djorgovski database jDaly & Djorgovsld 2004^ fitted to our model. The 
solid curve is the best fit curve with;^'^^^ =1.19 for 186 data points, and the best fit value is /fo = I.OO^qJJj, 



l |Daly & Djorgovski 2004^ developed a numerical method for a direct determination of the ex- 
pansion and acceleration rates, H{z) and q{z), from the data, just using the dimensionless coor- 
dinate distance y(z), without making any assumptions about the nature or evolution of the dark 
energy. They actually use the equation 



- q(z) = aalc? = 1 + (1 + z) {dy/dzr\d^y/dz^) , 



(43) 



valid for k - Q. Equation ( I43> depends only upon the Friedman-Robertson-Walker line element 
and the relation {I + z) - ao/a(t). Thus, this expression for q(z) is valid for any homogeneous 
and isotropic universe in which (1 + z) = a()/a{t), and is therefore quite general and can be 
compared with any model to account for the acceleration of the universe. This new approach 
has the advantage of being model independent, but it introduces larger errors in the estimation 
of q(z), because the numerical derivation is very sensitive to the size and quality of the data. An 
additional problem is posed by the sparse and not complete coverage of the z-range of interest. 
Measurement errors are propagated in the standard way leading to estimated uncertainties of the 
fitted values. In Fig. (I14> we compare the q{z) obtained by Daly & Djorgovski from their full 
data set with our best fit model. 



3.2. The Sunyaev-Zeldovich/X-ray method 

In this section we discuss how the parameters of our model can be also constrained by the an- 
gular diameter distance Da as measured using the Sunyaev-Zeldovich effect (SZE) and the ther- 
mal bremsstrahlung (X-ray brightness data) for galaxy clusters. Actually in a homogenous and 
isotropic cosmological model the angular diameter distance can be easily related to the coordi- 
nate distance leading to 



Da 



Ho (I 



+ z)Jo mo 



d^. 



(44) 
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Fig. 14. A first look at the allowed region for q(z), obtained by Daly & Djorgovski from the full data 
set (shadow area). Approximated polynomial corresponding to a z-window Az = 0.6 is shown with the 
black thin solid line, with the black thick dashed lines are shown the approximated polynomial fitted to 
the smoothed data at ±lcr range, and corresponding to a z- window Az = 0.4. The blue solid line shows 
the deceleration function, q(z), for our model corresponding to the the best fit values /fo = l.OO^jjJJj, a\ = 
2.5;^Q2. It is compared to the corresponding function for the standard ACDM model with Q.^ = 0.3, Q.\ = 
0.7. 



Distance measurements using SZE and X-ray emission from the intracluster medium are based 
on the fact that these processes depend on different combinations of some parameters of the 
clusters (see lBirkinshaw 1999l and references therein). The SZE is a result of the inverse Compton 
scattering of the CMB photons on hot electrons of the intracluster gas. The number of photons is 
preserved, but photons gain energy and thus a decrement of the temperature is generated in the 
Rayleigh-Jeans part of the black-body spectrum while an increment appears in the Wien region. 
We limit our analysis to the so called thermal or static SZE, which is present in all the clusters, 
neglecting the kinematic effect, which is present only in clusters with a nonzero peculiar velocity 
with respect to the Hubble flow along the line of sight. Typically the thermal SZE is an order of 
magnitude larger than the kinematic one. The shift of temperature is; 



To 



y 



xcoth 

hv 



©■ 



(45) 



where x - is a dimensionless variable, T is the shifted radiation temperature, Tq is the 

ksT 

unperturbed CMB temperature and y is the so called Compton parameter, defined as the optical 
depth T = (Tt- J ricdl times the energy gain per scattering: 

' kBTe 



y-' 



/: 



rrieCTTdl. 



(46) 



In Eq. ( I46> . T^ is the temperature of the electrons in the intracluster gas, me is the electron mass. 
He is the number density of the electrons, and ctt is the cross section of Thompson electron 
scattering. We have used the condition T,, » Tq (T^ is of the order of 10^ K and Tq, is the CMB 
temperature ^ 2.7K) and we assumed that the CMB temperature varies linearly with redshift 
what implies that after recombination the CMB radiation cools adiabaticaUy with no injection of 
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energy in the form of photons. In the low frequency regime of the Ray leigh- Jeans approximation 
we obtain 

^--2y. (47) 

The next step to quantify the SZE decrement is to specify the model for the intracluster electron 
density and temperature distribution. The most commonly used model is the so called isothermal 
jS model of lCavaliere & Fusco Femiano I In this model 



rieir) = n^J 1 + ( — 



2\-- 



(48) 
T,(r) = r,„ , (49) 



where «e„ and T,,^ are respectively the central electron number density and temperature of the 
intracluster electron gas, r^ andyS are fitting parameters connected with the model JSarazin 1988> . 
The relative temperature shift is given by 

AT IkRO-TTer, tier, 

1 nieC'- 



where 

dl, (51) 



f 

Jo 






which depends only on the geometry and the extension of the cluster along the line of sight. In 
Eq.j51>. I is the coordinate along the line of sight, r^ - f + R^, and R^ = x^ + y^. A simple geo- 
metrical argument converts the integral in Eg. (15 111 into an angular form. Introducing the angular 
diameter distance, cIa, to the cluster we can rewrite ( 15 Ot as 



AT(6 ^ 0) (TTkBTen^ ^r(l) c 

= -2 vn — -— — — di, (52) 

To me r(f)^o 

where T,, is the gas temperature. The factor — dA in Eq. ( I52t carries the dependence of the ther- 
mal SZE on the cosmological models (for a discussion of the dependence of dA on the standard 
ACDM model see dPemianski & al. 2003» . From Eq. ( I52K we also note that the central electron 
number density is proportional to the inverse of the angular diameter distance, actually 

<-^h "« 

^From an independent point of view the central number density of electrons can be also measured 
by fitting the X-ray surface brightness profile, Sx'^ J nlA(Te)dl, where the integration is along 
the line of sight and A(re) is the X-ray emissivity at the electron temperature T^. It turns out that 

«fo - JP. (54) 

V "A 

By ehminating WeO from Eqs. ( I54t . and ( I53> . one can solve for the angular diameter distance, 
yielding 

dA ^ — ■ (55) 
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Fig. 15. Observational SZE data fitted to our model with the best fit values Hq = 0.97;^j;j53, a\ 
and h = 0.75 ± 0.05 
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Recently distances to 18 clusters with redshift ranging from z ~ 0.14 to z ~ 0.78 have been 
determined from a likelihood joint analysis of SZE and X-ray observations (see Table 7 in 
IReese & al. 2002> . We perform our analysis using angular diameter distance measurements for a 
sample of 44 clusters, containing the 18 above mentioned clusters and other 24 known previously 
(see lBirkinshaw 1999I I. We perform a statistical analysis on the SZE data defining the following 
merit function : 



x\auHo) 



M-3 ^ 



{DA(zr,auHQ)-Dd 



o-i 



(56) 



We obtain ^2^^ =1.14 for 44 data points, and the best fit values are Hq = 0.97+°°^, ai = 3.2+°[. 
We also get h = 0.75 + 0.05. In Fig. ( I15> we compare the best fit curve with the observational 
SZE data. 



3.3. Gamma-Ray Burst Hubble Diagram 

Gamma-ray bursts (GRBs) are bright explosions visible across most of the Universe, certainly out 
to redshifts of z = 4.5 and likely out to z ~ 10. Recent studies have pointed out that GRBs may 
be used as standard cosmological candles. The prompt energy released during burst spans nearly 
three orders of magnitude, and the distribution of the opening angles of the emission, as deduced 
from the timing of the achromatic steepening of the afterglow emission, spans a similar wide 
range of values. However, when the apparently isotropic energy release and the conic opening 
of the emission are combined to infer the intrinsic, true energy release, the resulting distribution 
does not widen, as is expected for uncorrected data, but shrinks to a very well determined value 
JFrail & Kulkarni 2003> . with a remarkably small (one-sided) scattering, corresponding to about 
a factor of 2 in total energy. Similar studies in the X-ray band have reproduced the same results. It 
is thus very tempting to study to what extent this property of GRBs makes them suitable cosmo- 
logical standard candles. Schaefer JSchaefer 2003> proposed using two well known correlations 
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of the GRBs luminosity (with variabiHty, and with time delay) to the same end, while other ex- 
ploited the recently reported relationship between the beaming-corrected y-ray energy and the 
locally observed peak energy of GRBs (see for instance FPai & al. 2004> . As for the possible vari- 
ation of ambient density from burst to burst, which may widen the distribution of bursts energies. 
Frail & Kulkarni JFrail & Kulkarni 2003> remarked that this spread is already contained in their 
data sample, and yet the distribution of energy released is still very narrow. There are at least two 
reasons why GRBs are better than type la supernovae as cosmological candles. On the one hand, 
GRBs are easy to find and locate: even 1980s technology allowed BATSE to locate ~1 GRB per 
day, despite an incompleteness of about 1/3, making the build-up of a 300-object database a 
one-year enterprise. The Swift satellite launched on 20 November 2004, is expected to detect 
GRBs at about the same rate as BATSE, but with a nearly perfect capacity for identifying their 
redshifts simultaneously with the afterglow observations ^ . Second, GRBs have been detected 
out to very high redshifts: even the current sample of about 40 objects contains several events 
with z > 3, with one (GRB 000131) at z = 4.5. This should be contrasted with the difficulty of 
locating SN at z > 1, and the absolute lack of any SN with z > 2. On the other hand, the dis- 
tribution of luminosities of SNIa is narrower than the distribution of energy released by GRBs, 
corresponding to a magnitude dispersion ctm - 0.18 rather than ctm - 0.75. Thus GRBs may 
provide a complementary standard candle, out to distances which cannot be probed by SNIa, their 
major limitation being the larger intrinsic scatter of the energy released, as compared to the small 
scatter in peak luminosities of SNIa. There currently exists enough information to calibrate lu- 
minosity distances and independent redshifts for nine bursts (ISchaefer 2003> . These bursts were 
all detected by BATSE with redshifts measured from optical spectra of either the afterglow or 
the host galaxy. The highly unusual GRB980425 (associated with supernova SN1998bw) is not 
included because it is likely to be qualitatively diff'erent from the classical GRBs. Bursts with 
red shifts that were not recorded by BATSE cannot yet have their observed parameters converted 
to energies and fluxes that are comparable with BATSE data. We perform our analysis using the 
data shown in Fig. (I16> with the distance modulus ji, given by Eq. (I39> . To this aim, the only 
difference with respect to the SNIa is that we slightly modify the correction term of Eq. J38> . into 

AGeffit) 

^C'ff = 2.5y ^^y . (57) 

(InlO)Gf// 

We expect that y is of order unity, so that the G-correction would be roughly half a magnitude. 

We obiamxl^j = 1.09, and the best fit value is Hq = 1+°°^^ ay = 2.8+°'^, and h = 0.66 + 0.05, 

which are compatible with the SNIa results. We also confirm that y - 1.5 as in Eq. (I38> . In Fig. 

Ml\ we compare the best fit curve with both the GRBs and the SNIa Gold Sample. 

3.4. The gas fraction in clusters 

In this section we consider a recently proposed test based on the gas mass fraction in galaxy 
clusters ( JAUen & al. 2002> . Both theoretical arguments and numerical simulations predict that 

^ http://swift.gsfc.nasa.gov/docs/swift/proposals/appendix_f.html 
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Fig. 16. Hubble diagram for the BATSE gamma ray bursts JSchaefer 20031 up to z = 4.5. 
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Fig. 17. Observational Hubble diagram for the 
SNIa Gold Sample (Riess et al. 2004), the SNLS 
data JAstier & al. 2005 1 (filled boxes), and the BATSE 
GRBs data JSchaefer 2003> (empty boxes) fitted to our 
model. The solid curve is the best fit curve with Hq = 
l!ooI. f^i = 2.8;;;-^, and h = 0.68 ± 0.05. 




Fig. 18. The same as in the Eq. <17t , but zooming on the high 
redshift GRBs. 



the baryonic mass fraction in the largest relaxed galaxy clusters should not depend on the 
redshift, and should provide an estimate of the cosmological baryonic density parameter Q/, 
<Eke & al. 1998> . The baryonic content in galaxy clusters is dominated by the hot X-ray emit- 
ting intra-cluster gas so that what is actually measured is the gas mass fraction fgas and it is 
this quantity that should not depend on the redshift. Moreover, it is expected that the baryonic 
mass fraction in clusters equals the universal ratio Q.b/Q.M so that fgas should indeed be given 
by by.{D.bl^M), where Q^ is the matter density parameter, and the multiplicative factor b is 
motivated by simulations that suggest that the gas fraction is lower than the universal ratio be- 
cause of processes that convert part of the gas into stars or eject it out of the cluster altogether. 
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Following the procedure described in JAllen & al. 20021 IMen & al. 2004t . we adopt the SCDM 
model (i.e., a flat universe with Q^f = 1 and h - 0.5, where h is the Hubble constant in units of 
100 km s ' Mpc ') as a reference cosmology in making the measurements so that the theoretical 
expectation for the apparent variation of fgas with the redshift is: 



fgas(z) = 



(I +0.19 ^/h)£l 



M 



Z)f«^(z)l'-' 



'\z) 



(58) 



where D^j^^^ and Z)™°'' is the angular diameter distance for the SCDM and our model re- 
spectively. Allen & al. JAllen & al. 2002> have extensively analyzed the set of simulations 
in JEke & al. 1998^ to get b = 0.824±0.089, so in our analysis below, we set b = 0.824. Actually, 
we have checked that, for values in the 2cr range quoted above, the main results do not depend 
on b, as also on a^ . Moreover we have defined the following merit function : 

2_ 2 ( Q.bh^-0.02U \^ I h -0.72 \^ i b -0.824 \^ 
X -Xgas + y QQQ20 j ^\ 0.08 j ^\ 0.089 j ' ^ ' 

where we substitute the appropriate expression of ^m for our model 



1=1 



fgaAzi,auHo)-f^!;:(z,) 



cr„ 



(60) 



Here fgas(Zi) is the measured gas fraction in galaxy clusters at redshift z, with an error o-gi and 
the sum is over the Ngas clusters considered. Before presenting results of our analysis it is worth 
noting that, in order to estimate the gas fraction, it is necessary to evaluate the total cluster mass 
given by M,oi = -r^. Generally the standard assumption used to derive clusters masses from 
X-ray data is that the system is in hydrostatic equilibrium. This allows one to obtain a mass 
estimator just through the gas dynamical equilibrium equation: 



M(< r) 



Gj-mip 



dpgc 



d\nr 



(61) 



where ks is the Boltzmann constant, T the cluster gas temperature, // the mean molecular weight, 
Mp the proton mass, and pgas the gas mass density profile. Let us note that recently Allen & al. 
J Allen & al. 2004t have released a catalog of 26 large relaxed clusters with a precise measure- 
ment of both the gas mass fraction fgas and the redshift z- Actually to avoid possible systematic 
errors in the fgas measurement, it is desirable that the cluster is both highly luminous (so that 
the S/N ratio is high) and relaxed, so that both merging processes and cooling flows are absent. 
We use these data to perform our likelihood analysis, getting ;^'^ = 1.17 for 26 data points, and 
ai = 2.5!°-4^ //o = 0.98 ± 0.04, h = 0.72 ± 0.05, and w^ = -0.82 + 0.1. To complete our 
analysis we carry out a brief comparison of our results with similar recent results of Lima et al. 
( ILima & al. 2003t . where the equation of state characterizing the dark energy component is con- 
strained by using galaxy cluster x-ray data. In their analysis, however, they consider quintessence 
models in standard gravity theories, with a non evolving equation of state, but they allow the 
so-called phantom dark energy with w < -1, what violates the null energy condition. As best 
fit value of w to the data of ( lAUen & al. 2002> they obtain w - -1.29+Q^g2- In order to directly 
compare this result with our analysis we first fit the model considered in (ILima & al. 2003> to 
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Fig. 19. In the diagram we plot the best fit curve to the fgas data for our nmc model (green thin line) and 
for the quintessence model (black thick line) considered in ILima & al. 2003> . It is interesting to note the 
different behaviour of the two curves, even if the statistical significance of the best fit procedure for these 
two models is comparable: the best fit relative to our nmc model seems to be dominated by smaller redshift 
data, while the one relative to the Lima & al. model by higher redshift data. 



the updated and wider dataset of J Allen & al. 2004t . used in our analysis. To this aim we also 
refer to the model function fgasiz), and the merit function ;t'^, defined in the Eqs. J58l andl59t 
respectively. We get;^^ = 1.175 for 26 data points, and ^m = 0.23+|^°^, h = 0.76+°°^^ and 
w - -1.11 ± 0.35, so w < -1, what corresponds to a phantom energy. We note that our model, 
instead, gives w^ = -0.82 ±0.1, what does not violate the null energy condition. In Fig. ( I19> we 
compare the best fit curves for our and the Lima & al. model with the observational data. 



4. Growth of density perturbations 

In this section we consider the behavior of scalar density perturbations in the longitudinal gauge 
ds-^ = -(1 + 2<b)dt^ + a^(l - 2<t>)dx^. Actually, while in the framework of the minimally cou- 
pled theory we have to deal with a fully relativistic component, which becomes homogeneous 
on scales smaller than the horizon, so that standard quintessence cannot cluster on such scales. 
In the non minimally coupled quintessence theories it is possible to separate a pure gravita- 
tional term both in the stress-energy tensor Tf^y, and in the energy density p^ , so the situation 
changes, and it is necessary to consider also fluctuations of the scalar field. However, it turns 
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out (IBoisseau & al. 2000 I IRiazuelo & Uzan 2002> that the equation for dustlike matter density 
perturbations inside the horizon can be written as follows: 

Sm + 2H5m - -Gcav Pm 5„, ^ , (62) 

with Gcav is the effective gravitational constant between two test masses and is defined by 

_ 1 ( 2F + 4idF/d<p)^ \ 

The equation \62\ describes in the non minimally coupled models evolution of the CDM density 
contrast, 5,„ = 6p,„lpm, for perturbations inside the horizon. In our model the Ea.( l62t is rather 
complicated and takes the form 

3 (aj, + e"i' + a2t - 1) 

e"' [axa2 + e"' [a] - \6V(^) 

^ :;5m - 0, 

6 (q;2 + e"i + ffj - 1) (ffs + e"'' + 0-2? - 1) flf 

where 

ai(-a3+e"'+l)-3Ho(ai+e"' -1) 

Q'2 = ^^ ■ (65) 

3//0 + Qfi - 2 



The Eq. ( I64> does not admit exact solutions, and can be solved only numerically. However, since 
with our choice of normaUzation the whole history of the Universe is confined to the range 
f e [0, 1] and therefore to study the behavior of the solution for f ^ we can expand the expo- 
nential functions in Eq. \6A\ in series around f = 0. We obtain an integrable Fuchsian differential 
equation, which is a hypergeometric equation. We then use the obtained exact solution to set the 
initial conditions at f = to numerically integrate Eq. ( I64> in the whole range [0, 1]. We use the 
growing mode 6+ to construct the growth index / as 

/^^, (66) 

din a 

where a is the scale factor. Once we know how the growth index / evolves with redshift and 

how it depends on our model parameters, we can use the available observational data to estimate 

the values of such parameters, and the present value of Qmo- The 2dFGRS team has recently 

collected positions and redshifts of about 220000 galaxies and presented a detailed analysis of 

/ 
the two-point correlation function. They measured the redshift distortion parameter/? - — , where 

b 

b is the bias parameter describing the difference in the distribution of galaxies and mass, and 
obtained that yS|-.^o. 15 = 0.49 ± 0.09 and b - 1 .04 ± 0. 11 . From the observationally determined p 
and b it is now straightforward to get the value of the growth index at z = 0.15 corresponding to 
the effective depth of the survey. Verde & al. (IVerde & al. 200 U used the bispectrum of 2dFGRS 
galaxies, and Lahav & al. JLahav & al. 2002t combined the 2dFGRS data with CMB data, and 
they obtained 

b„erde = 1.04+0.11, (67) 

bu^Hav = 1.19 + 0.09. (68) 
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Fig. 20. The growth index / in different cosmological models: the solid line corresponds to our non mini- 
mally coupled model. The dashed curves correspond to a more standard quintessence model with an expo- 
nential potential (described in lDemianski & al. 2005> . 



Using these two values for b we calculated the value of the growth index / at z = 0.15, we get 
respectively 



/i = 0.51+0.1, 
/2 = 0.58 ±0.11. 



(69) 
(70) 



To evaluate the growth index at z = 0.15 we first have to invert the z - t relation and find fo.is- 
Then, substituting z = 0.15 and the two values of /i and /2 we calculate Hq and ai. Actually 
the z - t relation is rather involved and cannot be exactly inverted, so we apply this procedure 
numerically. We get ai = 3.3 + 0.05, Hq = 0.98+°°^, Vo = 0.5 + 0.06 which coiTesponds to 
D.Ao = 0.65 ± 0.08. In Fig. (I20> we show how the growth index is changing with redshift in 
our non minimally coupled model as compared with a standard quintessence model namely the 
minimally coupled exponential model described in JDemianski & al. 2005> . We note that at low 
redshift theoretical predictions of these different models are not distinguishable, independent 
measurements from large redshift surveys at different depths can disentangle this degeneracy. 



5. Conclusions 

In this paper we have shown that in the framework of non minimally coupled scalar tensor the- 
ory of gravitation it is possible to consider homogeneous and isotropic cosmological models 
with time dependent dark energy component. These models have a very interesting feature of 
producing in a natural way an epoch of accelerated expansion. In these models initially the over 
all density of the universe is dominated by matter (for the sake of simplicity we do not include 
radiation into our consideration) and later on the energy density of the scalar field becomes dom- 
inant and the universe enters an accelerated phase of its evolution. It turns out that at large t 
asymptotically w^ -^ -1, and such a transition to the asymptotical value is responsible for the 
accelerated expansion. The equation of state can admit a superquintessence behavior {w < -1), 
without violating the weak energy condition, just as an effect of the transition toward w ^ -1. 
The energy density of the scalar field decreases with time but before the present epoch it starts 
to dominate the expansion rate of the universe and asymptotically for f — > oo the universe is 
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reaching a de Sitter stage. To check viability of our model we have compared its predictions with 
the available observational data and it turned out that with an appropriate choice of parameters 
our model is reproducing the observed characteristics of the universe. Let us note that since in 
order to exactly solve the dynamical equations we have not included radiation into our consid- 
eration we do not use CMB data and observed abundances of Ught elements in confrontation of 
theoretical predictions of our model with observational data. 

In Tab.Q we present results of our analysis, they show that predictions of our model are fully 
compatible with the recent observational data. 

Table 1. The basic cosmological parameters derived from our model are compared with obser- 
vational data. 



o-i 


Ho 


"Ao 


w 


dataset 


2 Q+0.3 


1 Q+0.03 
^•'^-0.04 


a 04+O.O6 
"•°^-0.07 


-0.86 ± 0.06 


Gold SNIa + SNLS 


T 2+01 


A Q7+0.04 


0.62 ± 0.08 


-1.1 ±0.15 


SZe 


^•-'-0.2 


1 +0.05 
^-0.03 


0.8 ±0.15 


-0.76 ±0.1 


dimensionless coordinate 


r, O+O.l 


1 +0.05 
^-0.04 


0.84 ± 0.05 


-0.81 ± 0.07 


GRBs 


^•■'-0.1 


0.98 + 0.04 


0.76 ± 0.06 


-0.82 ±0.1 


gas fraction in clusters 


3.3 ± 0.05 


0.98!°:°^ 


0.65 ± 0.08 


-1.1 ±0.07 


galaxies peculiar velocity 


3.2 ± 0.04 


0.99 ± +0.02 


0.77 + 0.03 


-0.9 ± 0.04 


averaged mean 



Acknowledgments 

This work was supported in part by the grant of Polish Ministry of Science and Information 
Society Technologies 1-P03D-014-26, INFN Nal2 and the PRIN DRACO. The authors are very 
grateful to prof. Djorgovski for providing the data that we used in Sec. 3. 1 . 1 . Of course we take 
the full responsibility of the fitting procedure. 



6. Appendix I 

As was already noted the solution of the coupled system of the Einstein equations and generalized 

Klein-Gordon equation describing our model can be written in the form 

,?+] 



. „ , 2.t^+6.t+3 

a(t) = A{s) I B(s)t'+^ "'' 5^ I ' """ ' 



I Vq 3 D\ S" (2«3)2 

\ 7(s) lo/ 

where A(s), B(s), C(s), y(s) and;if(s) are given by 

,£i/(^ + 3)£\^ 



A(s) = C^C..))- 



3r(^) 



(71) 



(72) 



(73) 
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,_e+2s) 

(2..+3) (S + 3)I.\ 2<'+'* 

and 

2s + 3 
12(s + l){s + 2) 

^(*> = "T^ • (77> 

2s + 3 
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